Tracking SARS-CoV-2 introductions in Mozambique using pandemic-scale phylogenies: a retrospective observational study.

BACKGROUND: From the start of the SARS-CoV-2 outbreak, global sequencing efforts have generated an unprecedented amount of genomic data. Nonetheless, unequal sampling between high-income and low-income countries hinders the implementation of genomic surveillance systems at the global and local level. Filling the knowledge gaps of genomic information and understanding pandemic dynamics in low-income countries is essential for public health decision making and to prepare for future pandemics. In this context, we aimed to discover the timing and origin of SARS-CoV-2 variant introductions in Mozambique, taking advantage of pandemic-scale phylogenies. METHODS: We did a retrospective, observational study in southern Mozambique. Patients from Manhiça presenting with respiratory symptoms were recruited, and those enrolled in clinical trials were excluded. Data were included from three sources: (1) a prospective hospital-based surveillance study (MozCOVID), recruiting patients living in Manhiça, attending the Manhiça district hospital, and fulfilling the criteria of suspected COVID-19 case according to WHO; (2) symptomatic and asymptomatic individuals with SARS-CoV-2 infection recruited by the National Surveillance system; and (3) sequences from SARS-CoV-2-infected Mozambican cases deposited on the Global Initiative on Sharing Avian Influenza Data database. Positive samples amenable for sequencing were analysed. We used Ultrafast Sample placement on Existing tRees to understand the dynamics of beta and delta waves, using available genomic data. This tool can reconstruct a phylogeny with millions of sequences by efficient sample placement in a tree. We reconstructed a phylogeny (~7·6 million sequences) adding new and publicly available beta and delta sequences. FINDINGS: A total of 5793 patients were recruited between Nov 1, 2020, and Aug 31, 2021. During this time, 133 328 COVID-19 cases were reported in Mozambique. 280 good quality new SARS-CoV-2 sequences were obtained after the inclusion criteria were applied and an additional 652 beta (B.1.351) and delta (B.1.617.2) public sequences were included from Mozambique. We evaluated 373 beta and 559 delta sequences. We identified 187 beta introductions (including 295 sequences), divided in 42 transmission groups and 145 unique introductions, mostly from South Africa, between August, 2020 and July, 2021. For delta, we identified 220 introductions (including 494 sequences), with 49 transmission groups and 171 unique introductions, mostly from the UK, India, and South Africa, between April and November, 2021. INTERPRETATION: The timing and origin of introductions suggests that movement restrictions effectively avoided introductions from non-African countries, but not from surrounding countries. Our results raise questions about the imbalance between the consequences of restrictions and health benefits. This new understanding of pandemic dynamics in Mozambique can be used to inform public health interventions to control the spread of new variants. FUNDING: European and Developing Countries Clinical Trials, European Research Council, Bill & Melinda Gates Foundation, and Agència de Gestió d'Ajuts Universitaris i de Recerca.

Deciphering the Tangible Spatio-Temporal Spread of a 25-Year Tuberculosis Outbreak Boosted by Social Determinants.

Outbreak strains of Mycobacterium tuberculosis are promising candidates as targets in the search for intrinsic determinants of transmissibility, as they are responsible for many cases with sustained transmission; however, the use of low-resolution typing methods and restricted geographical investigations represent flaws in assessing the success of long-lived outbreak strains. We can now address the nature of outbreak strains by combining large genomic data sets and phylodynamic approaches. We retrospectively sequenced the whole genome of representative samples assigned to an outbreak circulating in the Canary Islands (the GC strain) since 1993, which accounts for ~20% of local tuberculosis cases. We selected a panel of specific single nucleotide polymorphism (SNP) markers for an in-silico search for additional outbreak-related sequences within publicly available tuberculosis genomic data. Using this information, we inferred the origin, spread, and epidemiological parameters of the GC strain. Our approach allowed us to accurately trace the historical and more recent dispersion of the GC strain. We provide evidence of a highly successful nature within the Canarian archipelago but limited expansion abroad. Estimation of epidemiological parameters from genomic data disagree with a distinctive biology of the GC strain. With the increasing availability of genomic data allowing for the accurate inference of strain spread and critical epidemiological parameters, we can now revisit the link between Mycobacterium tuberculosis genotypes and transmission, as is routinely carried out for SARS-CoV-2 variants of concern. We demonstrate that social determinants rather than intrinsically higher bacterial transmissibility better explain the success of the GC strain. Importantly, our approach can be used to trace and characterize strains of interest worldwide. IMPORTANCE Infectious disease outbreaks represent a significant problem for public health. Tracing outbreak expansion and understanding the main factors behind emergence and persistence remain critical to effective disease control. Our study allows researchers and public health authorities to use Whole-Genome Sequencing-based methods to trace outbreaks, and shows how available epidemiological information helps to evaluate the factors underpinning outbreak persistence. Taking advantage of all the freely available information placed in public repositories, researchers can accurately establish the expansion of an outbreak beyond original boundaries, and determine the potential risk of a strain to inform health authorities which, in turn, can define target strategies to mitigate expansion and persistence. Finally, we show the need to evaluate strain transmissibility in different geographic contexts to unequivocally associate spread to local or pathogenic factors, an important lesson taken from genomic surveillance of SARS-CoV-2.

Cumulative incidence of SARS-CoV-2 infection in the general population of the Valencian Community (Spain) after the surge of the Omicron BA.1 variant.

Studies investigating the cumulative incidence of and immune status against SARS-CoV-2 infection provide valuable information for shaping public health decision-making. A cross-sectional study on 935 participants, conducted in the Valencian Community (VC), measuring anti-SARS-CoV-2-receptor binding domain-RBD-total antibodies and anti-Nucleocapsid (N)-IgGs via electrochemiluminescence assays. Quantitation of neutralizing antibodies (NtAb) against ancestral and Omicron BA.1 and BA.2 variants and enumeration of SARS-CoV-2-S specific-IFNγ-producing CD4+ and CD8+ T cells was performed in 100 and 137 participants, respectively. The weighted cumulative incidence was 51.9% (95% confidence interval [CI]: 48.7-55.1) and was inversely related to age. Anti-RBD total antibodies were detected in 97% of participants; vaccinated and SARS-CoV-2-experienced (VAC-ex; n = 442) presented higher levels (p < 0.001) than vaccinated/naïve (VAC-n; n = 472) and nonvaccinated/experienced (UNVAC-ex; n = 63) subjects. Antibody levels correlated inversely with time elapsed since last vaccine dose in VAC-n (Rho, -0.52; p < 0.001) but not in VAC-ex (rho -0.02; p = 0.57). Heterologous booster shots resulted in increased anti-RBD antibody levels compared with homologous schedules in VAC-n, but not in VAC-ex. NtAbs against Omicron BA.1 were detected in 94%, 75%, and 50% of VAC-ex, VAC-n and UNVAC-ex groups, respectively. For Omicron BA.2, the figures were 97%, 84%, and 40%, respectively. SARS-CoV-2-S-reactive IFN-γ T cells were detected in 73%, 75%, and 64% of VAC-ex, VAC-n and UNVAC-ex, respectively. Median frequencies for both T-cell subsets were comparable across groups. In summary, by April 2022, around half of the VC population had been infected with SARS-CoV-2 and, due to extensive vaccination, displayed hybrid immunity.

Prevalence of SARS-CoV-2 and co-infection with malaria during the first wave of the pandemic (the Burkina Faso case).

Africa accounts for 1.5% of the global coronavirus disease 2019 (COVID-19) cases and 2.7% of deaths, but this low incidence has been partly attributed to the limited testing capacity in most countries. In addition, the population in many African countries is at high risk of infection with endemic infectious diseases such as malaria. Our aim is to determine the prevalence and circulation of SARS-CoV-2 variants, and the frequency of co-infection with the malaria parasite. We conducted serological tests and microscopy examinations on 998 volunteers of different ages and sexes in a random and stratified population sample in Burkina-Faso. In addition, nasopharyngeal samples were taken for RT-qPCR of SARS-CoV-2 and for whole viral genome sequencing. Our results show a 3.2 and a 2.5% of SARS-CoV-2 seroprevalence and PCR positivity; and 22% of malaria incidence, over the sampling period, with marked differences linked to age. Importantly, we found 8 cases of confirmed co-infection and 11 cases of suspected co-infection mostly in children and teenagers. Finally, we report the genome sequences of 13 SARS-CoV-2 isolates circulating in Burkina Faso at the time of analysis, assigned to lineages A.19, A.21, B.1.1.404, B.1.1.118, B.1 and grouped into clades; 19B, 20A, and 20B. This is the first population-based study about SARS-CoV-2 and malaria in Burkina Faso during the first wave of the pandemic, providing a relevant estimation of the real prevalence of SARS-CoV-2 and variants circulating in this Western African country. Besides, it highlights the non-negligible frequency of co-infection with malaria in African communities.

Evaluation of two different concentration methods for surveillance of human viruses in sewage and their effects on SARS-CoV-2 sequencing.

During the current COVID-19 pandemic, wastewater-based epidemiology (WBE) emerged as a reliable strategy both as a surveillance method and a way to provide an overview of the SARS-CoV-2 variants circulating among the population. Our objective was to compare two different concentration methods, a well-established aluminum-based procedure (AP) and the commercially available Maxwell® RSC Enviro Wastewater TNA Kit (TNA) for human enteric virus, viral indicators and SARS-CoV-2 surveillance. Additionally, both concentration methods were analyzed for their impact on viral infectivity, and nucleic acids obtained from each method were also evaluated by massive sequencing for SARS-CoV-2. The percentage of SARS-CoV-2 positive samples using the AP method accounted to 100 %, 83.3 %, and 33.3 % depending on the target region while 100 % positivity for these same three target regions was reported using the TNA procedure. The concentrations of norovirus GI, norovirus GII and HEV using the TNA method were significantly greater than for the AP method while no differences were reported for rotavirus, astrovirus, crAssphage and PMMoV. Furthermore, TNA kit in combination with the Artic v4 primer scheme yields the best SARS-CoV-2 sequencing results. Regarding impact on infectivity, the concentration method used by the TNA kit showed near-complete lysis of viruses. Our results suggest that although the performance of the TNA kit was higher than that of the aluminum procedure, both methods are suitable for the analysis of enveloped and non-enveloped viruses in wastewater by molecular methods.

Epidemiological and Genomic Analysis of a Large SARS-CoV-2 Outbreak in a Long-Term Care Facility in Catalonia, Spain.

Limiting outbreaks in long-term care facilities (LTCFs) is a cornerstone strategy to avoid an excess of COVID-19-related morbidity and mortality and to reduce its burden on the health system. We studied a large outbreak that occurred at an LTCF, combining methods of classical and genomic epidemiology analysis. The outbreak lasted for 31 days among residents, with an attack rate of 98% and 57% among residents and staff, respectively. The case fatality rate among residents was 16% (n = 15). Phylogenetic analysis of 59 SARS-CoV-2 isolates revealed the presence of two closely related viral variants in all cases (B.1.177 lineage), revealing a far more complex outbreak than initially thought and suggesting an initial spread driven by staff members. In turn, our results suggest that resident relocations to mitigate viral spread might have increased the risk of infection for staff members, creating secondary chains of transmission that were responsible for prolonging the outbreak. Our results highlight the importance of considering unnoticed chains of transmission early during an outbreak and making an adequate use and interpretation of diagnostic tests. Outbreak containment measures should be carefully tailored to each LTCF. IMPORTANCE The impact of COVID-19 on long-term care facilities (LTCFs) has been disproportionately large due to the high frailty of the residents. Here, we report epidemiological and genomic findings of a large outbreak that occurred at an LTCF, which ultimately affected almost all residents and nearly half of staff members. We found that the outbreak was initially driven by staff members; however, later resident relocation to limit the outbreak resulted in transmission from residents to staff members, evidencing the complexity and different phases of the outbreak. The phylogenetic analysis of SARS-CoV-2 isolates indicated that two closely related variants were responsible for the large outbreak. Our study highlights the importance of combining methods of classical and genomic epidemiology to take appropriate outbreak containment measures in LTCFs.

Correction: The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

[This corrects the article DOI: 10.1371/journal.ppat.1010631.].

SARS-CoV-2 Outbreak on a Spanish Mink Farm: Epidemiological, Molecular, and Pathological Studies.

Farmed minks have been reported to be highly susceptible to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and may represent a risk to humans. In this study, we describe the first outbreak of SARS-CoV-2 occurred on a mink farm in Spain, between June and July 2020, involving 92,700 animals. The outbreak started shortly after some farm workers became seropositive for SARS-CoV-2. Minks showed no clinical signs compatible with SARS-CoV-2 infection throughout the outbreak. Samples from 98 minks were collected for histopathological, serological, and molecular studies. Twenty out of 98 (20.4%) minks were positive by RT-qPCR and 82 out 92 (89%) seroconverted. This finding may reflect a rapid spread of the virus at the farm with most of the animals overcoming the infection. Additionally, SARS-CoV-2 was detected by RT-qPCR in 30% of brain samples from positive minks. Sequencing analysis showed that the mink sequences were not closely related with the other mink SARS-CoV-2 sequences available, and that this mink outbreak has its probable origin in one of the genetic variants that were prevalent in Spain during the first COVID-19 epidemic wave. Histological studies revealed bronchointerstitial pneumonia in some animals. Immunostaining of viral nucleocapsid was also observed in nasal turbinate tissue. Farmed minks could therefore constitute an important SARS-CoV-2 reservoir, contributing to virus spread among minks and humans. Consequently, continuous surveillance of mink farms is needed.

The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

The S:A222V point mutation, within the G clade, was characteristic of the 20E (EU1) SARS-CoV-2 variant identified in Spain in early summer 2020. This mutation has since reappeared in the Delta subvariant AY.4.2, raising questions about its specific effect on viral infection. We report combined serological, functional, structural and computational studies characterizing the impact of this mutation. Our results reveal that S:A222V promotes an increased RBD opening and slightly increases ACE2 binding as compared to the parent S:D614G clade. Finally, S:A222V does not reduce sera neutralization capacity, suggesting it does not affect vaccine effectiveness.

A Genomic Snapshot of the SARS-CoV-2 Pandemic in the Balearic Islands.

Objective: To analyze the SARS-CoV-2 genomic epidemiology in the Balearic Islands, a unique setting in which the course of the pandemic has been influenced by a complex interplay between insularity, severe social restrictions and tourism travels. Methods: Since the onset of the pandemic, more than 2,700 SARS-CoV-2 positive respiratory samples have been randomly selected and sequenced in the Balearic Islands. Genetic diversity of circulating variants was assessed by lineage assignment of consensus whole genome sequences with PANGOLIN and investigation of additional spike mutations. Results: Consensus sequences were assigned to 46 different PANGO lineages and 75% of genomes were classified within a VOC, VUI, or VUM variant according to the WHO definitions. Highest genetic diversity was documented in the island of Majorca (42 different lineages detected). Globally, lineages B.1.1.7 and B.1.617.2/AY.X were identified as the 2 major lineages circulating in the Balearic Islands during the pandemic, distantly followed by lineages B.1.177/B.1.177.X. However, in Ibiza/Formentera lineage distribution was slightly different and lineage B.1.221 was the third most prevalent. Temporal distribution analysis showed that B.1 and B.1.5 lineages dominated the first epidemic wave, lineage B.1.177 dominated the second and third, and lineage B.1.617.2 the fourth. Of note, lineage B.1.1.7 became the most prevalent circulating lineage during first half of 2021; however, it was not associated with an increased in COVID-19 cases likely due to severe social restrictions and limited travels. Additional spike mutations were rarely documented with the exception of mutation S:Q613H which has been detected in several genomes (n = 25) since July 2021. Conclusion: Virus evolution, mainly driven by the acquisition and selection of spike substitutions conferring biological advantages, social restrictions, and size population are apparently key factors for explaining the epidemic patterns registered in the Balearic Islands.

Severe Acute Respiratory Syndrome Coronavirus 2 Adaptive Immunity in Nursing Home Residents Following a Third Dose of the Comirnaty Coronavirus Disease 2019 Vaccine.

A third Comirnaty vaccine dose increased severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) receptor-binding domain antibody levels (median, 93-fold) and neutralizing antibody titers against Wuhan-Hu-1 (median, 57-fold), Beta (me 22-fold), Delta, (median, 43-fold), and Omicron (median, 8-fold) variants, but had less impact on S-reactive T-cell immunity in nursing home residents.

SARS-CoV-2 Omicron BA.1 variant breakthrough infections in nursing home residents after an homologous third dose of the Comirnaty® COVID-19 vaccine: Looking for correlates of protection.

We investigated whether peripheral blood levels of SARS-CoV-2 Spike (S) receptor binding domain antibodies (anti-RBD), neutralizing antibodies (NtAb) targeting Omicron S, and S-reactive-interferon (IFN)-γ-producing CD4+ and CD8+ T cells measured after a homologous booster dose (3D) with the Comirnaty® vaccine was associated with the likelihood of subsequent breakthrough infections due to the Omicron variant. An observational study including 146 nursing home residents (median age, 80 years; range, 66-99; 109 female) evaluated for an immunological response after 3D (at a median of 16 days). Anti-RBD total antibodies were measured by chemiluminescent immunoassay. NtAb were quantified by an Omicron S pseudotyped virus neutralization assay. SARS-CoV-2-S specific-IFNγ-producing CD4+ and CD8+ T cells were enumerated by whole-blood flow cytometry for intracellular cytokine staining. In total, 33/146 participants contracted breakthrough Omicron infection (symptomatic in 30/33) within 4 months after 3D. Anti-RBD antibody levels were comparable in infected and uninfected participants (21 123 vs. 24 723 BAU/ml; p = 0.34). Likewise, NtAb titers (reciprocal IC50 titer, 157 vs. 95; p = 0.32) and frequency of virus-reactive CD4+ (p = 0.82) and CD8+ (p = 0.91) T cells were similar across participants in both groups. anti-RBD antibody levels and NtAb titers estimated at around the time of infection were also comparable (3445 vs. 4345 BAU/ml; p = 0.59 and 188.5 vs. 88.9; p = 0.70, respectively). Having detectable NtAb against Omicron or SARS-CoV-2-S-reactive-IFNγ-producing CD4+ or CD8+ T cells after 3D was not correlated with increased protection from breakthrough infection (OR, 1.50; p = 0.54; OR, 0.0; p = 0.99 and OR 3.70; p = 0.23, respectively). None of the immune parameters evaluated herein, including NtAb titers against the Omicron variant, may reliably predict at the individual level the risk of contracting COVID-19 due to the Omicron variant in nursing home residents.

SARS-CoV-2 Delta-variant breakthrough infections in nursing home residents at midterm after Comirnaty® COVID-19 vaccination.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Delta variant breakthrough infections in nursing home residents following vaccination with Comirnaty® COVID-19 vaccine were characterized. In total, 201 participants (median age, 87 years; range, 64-100; 133 female) from two nursing homes in the Valencian community (Spain) were included. SARS-CoV-2-Spike (S) antibody responses were determined by a lateral flow immunocromatography (LFIC) assay and by quantitative electrochemiluminescent assay in LFIC-negative participants. SARS-CoV-2-S-IFNγ T cells were enumerated by flow cytometry in 10 participants. Nasopharyngeal SARS-CoV-2 RNA loads were quantified by real-time polymerase chain reaction assays. Vaccine breakthrough COVID-19 due to the Delta variant occurred in 39 residents (median age, 87 years; range, 69-96; 31 female) at a median of 6.5 months after vaccination (nine requiring hospitalization). Breakthrough infections occurred at a higher rate (p < 0.0001) in residents who had not been previously infected with SARS-CoV-2 (naïve) (33/108; 18%) than in those with prior diagnosis of SARS-CoV-2 infection (experienced) (6/93; 6.4%), and were more likely (p < 0.0001) to develop in residents who tested negative by LFIC (20/49) at 3 months after vaccination as compared to their LFIC-positive counterparts (19/142). Among LFIC-negative residents, a trend towards lower plasma anti-RBD antibody levels was noticed in those developing breakthrough infection (p = 0.16). SARS-CoV-2 RNA loads in nasopharyngeal specimens were lower in SARS-CoV-2-experienced residents (p < 0.001) and in those testing positive by LFIC (p = 0.13). The frequency of SARS-CoV-2-S-reactive T cells at 3 months was similar in LFIC-negative residents with (n = 7) or without (n = 3) breakthrough infection. Prior history of SARS-CoV-2 infection and detection of S-reactive antibodies by LFIC at 3 months is associated with a lower risk of Delta-variant breakthrough infection in nursing home residents at midterm after Comirnaty® COVID-19 vaccination.

SARS-CoV-2 adaptive immunity in nursing home residents up to eight months after two doses of the Comirnaty® COVID-19 vaccine.

Burgos JS (General Directorate of Research and Healthcare Supervision, Department of Health, Valencia Government, Valencia, Spain); Meneu de Guillerna R (Vice-President Foundation Research Institute in Public Services, Valencia, Spain); Vanaclocha Luna H (General Directorate of Public Health, Department of Health, Valencia Government, Valencia, Spain); Burks DJ (The Prince Felipe Research Center-CIPF-, Valencia, Spain; Cervantes A (INCLIVA Health Research Institute, Valencia, Spain); Comas I (Biomedicine Institute of Valencia, Spanish Research Council (CSIC); Díez-Domingo J (Foundation for the promotion of health and biomedical research of the Valencian Community-FISABIO-, Valencia, Spain); Peiro S (Foundation for the promotion of health and biomedical research of the Valencian Community-FISABIO-, Valencia, Spain); González-Candelas F (CIBER in Epidemiology and Public Health, Spain; Joint Research Unit "Infection and Public Health" FISABIO-University of Valencia, Valencia, Spain; Institute for Integrative Systems Biology (I2SysBio), CSIC-University of Valencia, Valencia, Spain); Ferrer Albiach C (Fundación Hospital Provincial de Castelló); Hernández-Aguado I (University Miguel Hernández, Alicante, Spain); Oliver Ramírez N (DataPop Alliance); Sánchez-Payá J (Preventive Medicine Service, Alicante General and University Hospital, Alicante, Spain; Alicante Institute of Health and Biomedical Research (ISABIAL), Alicante, Spain; Vento Torres M (Instituto de Investigación Sanitaria La Fe); Zapater Latorre E (Fundación Hospital General Universitario de València); Navarro D (Microbiology Service, Clinic University Hospital, INCLIVA Health Research Institute, Valencia, Spain;Department of Microbiology, School of Medicine, University of Valencia, Valencia, Spain).

Spatial and temporal distribution of SARS-CoV-2 diversity circulating in wastewater.

Wastewater-based epidemiology (WBE) has proven to be an effective tool for epidemiological surveillance of SARS-CoV-2 during the current COVID-19 pandemic. Furthermore, combining WBE together with high-throughput sequencing techniques can be useful for the analysis of SARS-CoV-2 viral diversity present in a given sample. The present study focuses on the genomic analysis of SARS-CoV-2 in 76 sewage samples collected during the three epidemiological waves that occurred in Spain from 14 wastewater treatment plants distributed throughout the country. The results obtained demonstrate that the metagenomic analysis of SARS-CoV-2 in wastewater allows the detection of mutations that define the B.1.1.7 lineage and the ability of the technique to anticipate the detection of certain mutations before they are detected in clinical samples. The study proves the usefulness of sewage sequencing to track Variants of Concern that can complement clinical testing to help in decision-making and in the analysis of the evolution of the pandemic.

Systematic Genomic and Clinical Analysis of Severe Acute Respiratory Syndrome Coronavirus 2 Reinfections and Recurrences Involving the Same Strain.

Estimates of the burden of severe acute respiratory syndrome coronavirus 2 reinfections are limited by the scarcity of population-level studies incorporating genomic support. We conducted a systematic study of reinfections in Madrid, Spain, supported by genomic viral analysis and host genetic analysis, to cleanse laboratory errors and to discriminate between reinfections and recurrences involving the same strain. Among the 41,195 cases diagnosed (March 2020-March 2021), 93 (0.23%) had 2 positive reverse transcription PCR tests (55-346 days apart). After eliminating cases with specimens not stored, of suboptimal sequence quality, or belonging to different persons, we obtained valid data from 22 cases. Of those, 4 (0.01%) cases were recurrences involving the same strain; case-patients were 39-93 years of age, and 3 were immunosuppressed. Eighteen (0.04%) cases were reinfections; patients were 19-84 years of age, and most had no relevant clinical history. The second episode was more severe in 8 cases.

Evolutionary and Phenotypic Characterization of Two Spike Mutations in European Lineage 20E of SARS-CoV-2.

We have detected two mutations in the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at amino acid positions 1163 and 1167 that appeared independently in multiple transmission clusters and different genetic backgrounds. Furthermore, both mutations appeared together in a cluster of 1,627 sequences belonging to clade 20E. This cluster is characterized by 12 additional single nucleotide polymorphisms but no deletions. The available structural information on the S protein in the pre- and postfusion conformations predicts that both mutations confer rigidity, which could potentially decrease viral fitness. Accordingly, we observed reduced infectivity of this spike genotype relative to the ancestral 20E sequence in vitro, and the levels of viral RNA in nasopharyngeal swabs were not significantly higher. Furthermore, the mutations did not impact thermal stability or antibody neutralization by sera from vaccinated individuals but moderately reduce neutralization by convalescent-phase sera from the early stages of the pandemic. Despite multiple successful appearances of the two spike mutations during the first year of SARS-CoV-2 evolution, the genotype with both mutations was displaced upon the expansion of the 20I (Alpha) variant. The midterm fate of the genotype investigated was consistent with the lack of advantage observed in the clinical and experimental data. IMPORTANCE We observed repeated, independent emergence of mutations in the SARS-CoV-2 spike involving amino acids 1163 and 1167, within the HR2 functional motif. Conclusions derived from evolutionary and genomic diversity analysis suggest that the co-occurrence of both mutations might pose an advantage for the virus and therefore a threat to effective control of the epidemic. However, biological characterization, including in vitro experiments and analysis of clinical data, indicated no clear benefit in terms of stability or infectivity. In agreement with this, continuous epidemiological surveillance conducted months after the first observations revealed that both mutations did not successfully outcompete other variants and stopped circulating 9 months after their initial detection. Additionally, we evaluated the potential of both mutations to escape neutralizing antibodies, finding that the presence of these two mutations on their own is not likely to confer antibody escape. Our results provide an example of how newly emerged spike mutations can be assessed to better understand the risk posed by new variants and indicate that some spike mutations confer no clear advantage to the virus despite independently emerging multiple times and are eventually displaced by fitter variants.